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Domain Wall Fermions and chiral symmetry restoration rate. 
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Domain Wall Fermions utilize an extra space time dimension to provide a method for restoring the regularization 
induced chiral symmetry breaking in lattice vector gauge theories even at finite lattice spacing. The breaking is 
restored at an exponential rate as the size of the extra dimension increases. As a precursor to lattice QCD studies 
the dependence of the restoration rate to the other parameters of the theory and, in particular, the lattice spacing 
is investigated in the context of the two flavor lattice Schwinger model. 



A massive vector theory in 2n + 1 dimensions 
with a mass term m(s) that depends on the peri- 
odic 2n + 1 dimension s such that m(s) = toq for 
1 < s < Ls/2 and m(s) = —mo for Ls/2 < s < 
Ls, where Lg is the size of the 2n + 1 dimension, 
develops a Dirac fermion in 2n dimensions with 
chiral components separated and localized along 
the 2n + 1 direction s at s = and s = Ls/2 [Q. 
Perturbative calculations indicate exponential lo- 
calization along s with the mixing of the chiral 
components vanishing exponentially fast with Ls. 
At Ls — cx) the mixing is zero and the regulariza- 
tion leaves intact the chiral symmetry. 

The Ls = CX) case can be handled using the 
Overlap formalism Q|. Therefore the Overlap is 
an ideal regulator for vector theories. However, 
calculations using the Overlap involve diagonal- 
ization and determinant evaluation of matrices 
with size ~ volume x volume. Because of this 
large scale dynamical QCD calculations with the 
Overlap are beyond the capacity of present day 
supercomputers (however, see |^ for a promising 
new development). 

An alternative is to keep Ls finite and view it 
as an extra parameter that controls the regular- 
ization induced chiral symmetry breaking. This 
formulation goes under the name of Domain Wall 
Fermions (DWF). Unlike Wilson or staggered 
fermions the chiral limit can be approached even 
at finite lattice spacing by increasing Ls- The 
theory can be simulated using standard Hybrid 
Monte Carlo techniques. However, before DWF 



can be used in dynamical simulations of QCD the 
following questions should be answered: 

1) Is the regularization induced chiral symme- 
try breaking restored exponentially fast as Ls in- 
creases in the full dynamical theory? 

2) How does the functional behavior of the 
restoration and the restoration rate depend on 
the lattice spacing a, the fermion mass m/ and 
the domain wall height toq? 

3) For the interesting a and uif how large should 
Ls be for the effects of the regularization induced 
chiral symmetry breaking to be negligible? 

This work addresses these questions in the con- 
text of the Nf = 2 flavor lattice Schwinger model 
at fixed physical volume. For a detailed analysis 
see For some interesting work that also par- 
tially addresses these questions see A variant 
domain wall model is used where m(s) = mo 
for I < s < Ls but the boundary conditions along 
the 2n-|-l direction are free. In this model the chi- 
ral components are localized at s = 1 and s — Ls 
and as a result the decay region is twice as long 
for the same Ls. Also, an explicit coupling of 
the two chiral components at s = and s — Ls 
with strength m/ is introduced to allow for linear 
control over the fermion mass. The action is the 
same as the one for a2n + l Wilson type fermion 
except that the gauge field "lives" only on the 2n 
dimensional lattice, the diagonal term is of the 
form 2n + 1 — mo and the boundary conditions 
along the 2n + 1 dimension are free. 

An explicit calculation B of the free theory 
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propagator gives the mass of the hghtest mode: 

"icff = m-o(2 - mo) [to/ + (1 - toq)"^"] 

This strongly suggests the pattern for chiral sym- 
metry restoration in the model. Recent calcula- 
tions of the one loop correction result in a renor- 
malization of toq . 

Two observables are used to probe the sym- 
metries of the theory: the chiral condensate 
(V'V') and the t'Hooft vertex (w) where w = 

YlfJi '^r'^l + riilfi i^Li^R- The chiral condensate 
is used to probe chiral symmetry. At toj = 
and Ls = oo (tpip) = 0- When Lg is finite 
the regularization breaks chiral symmetry and 
{tpt/j) ^ 0. From perturbative calculations one 
expects that (V'V') should approach zero exponen- 
tially fast with increasing Lg. The t' Hooft vertex 
is used to probe the anomalously broken U{1) ax- 
ial symmetry. At nif = 0, Lg = oo one should 
find that (w) 0. 

In order to study the approach to the Lg ^ oo 
limit it is very useful to calculate the values of 
these observables at = oo using the Overlap 
formalism. In two dimensions this is not compu- 
tationally demanding and can be done as in |^]. 
The results, interesting in their own right, are 
shown in figures 1 and 2. In figure 1 {ipip)/my is 
plotted vs. To/ for L — 6, toq — 0.9 and = 
3.0. go is the coupling constant ^il = g^L/^, 
m.y — ^y/Nf is the photon mass, L is the lattice 
size in lattice units and / in physical units. The 
fits are to (ipip) = AmF^. Both fits have a per 
degree of freedom of about one. For to/ < 0.1 
p = 0.996(3), while for to/ > 0.1 p = 0.32(2). 
This is in agreement with the analytical solutions 
found in [|lO|. In figure 2 < w > /to^ is plotted 
vs. TO/ for the same parameters. The dotted lines 
are the to/ = result ± the statistical error. As 
expected at to/ = 0, {w) ^ 0. 

The finite Lg model is studied using a stan- 
dard Hybrid Monte Carlo algorithm. In order to 
reveal the chiral symmetry breaking due to the 
DWF regularization, {ipip) is calculated for var- 
ious Lg at TO/ =0. It is plotted in figure 3 for 
fixed physical volume /iZ = 3.0, toq = 0.9 and four 
lattice spacings fjd/L = /la, with L = 6,8, 10, 12 
corresponding to the lines from top to bottom. 
The data is consistent with exponential decay in- 



dicated by the fitted lines (x^ per degree of free- 
dom « 1). However, at the larger spacing L — Q 
the data also fits well to a power law behavior. 
This is not so at the smallest spacing L = 12 (a 
power law fit has per degree of freedom « 31). 




Figure 1. Chiral condensate vs. fermion mass. 
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Figure 2. t' Hooft vertex vs. fermion mass. 




Figure 3. Chiral condensate vs. Lg for to/ = 0. 

This behavior is consistent with a picture where 
chiral symmetry is restored with a fast exponen- 
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tial decay rate for up to some value and with 
a slower exponential decay rate for Lg above that 
value. For the range of lattice spacings used, the 
inflection appeared at Lg « 10. Using a sim- 
ple model, it was found in [|| that the first fast 
decay can be associated with restoration of chi- 
ral symmetry in the zero topological sector while 
the second slower decay can be associated with 
the regions of gauge field configuration space that 
connect the q = and q — ±.l topological sectors. 
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Figure 4. Chiral condensate vs. Lg ioi irif ^ 0. 
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Figure 5. t' Hooft vertex vs. Lg for my 7^ 0. 

The effects of the size of the lattice spacing a 
to the two decays are apparent. The fast decay 
becomes faster as a is decreased. The vanishing 
of {tpip) is consistent with a form e^'^^" with 
being roughly a linear function of a. However, 
more data at smaller a are needed before one can 
be confident that this is the correct scaling form. 
The second slower decay also becomes faster as 
a is decreased and it differs less from the faster 
decay as a becomes smaller. This behavior can 
also be understood using the simple model in . 



For small but non zero nif the values of ('0V') 
and (w) are presented in figures 4 and 5 for two 
different lattice spacings set by L = 4 (lower 
curve) and L = 10 (upper curve) at toq = 0.9. 
The physical volume and m/L are fixed at fil = 
3.0 and m/i = 2.0. The fits are to A + Be'"^' . 
The dotted lines are the Lg = 00 results ± the er- 
ror. The cross is the coefficient A. The Lg ~ 00 
numbers are approached in a way that is con- 
sistent with exponential decay with a rate that 
becomes faster as the lattice spacing decreases. 
Finally, it was also found in Q that the larger the 
fermion mass the sooner the Lg — 00 value was 
approached (within a few percent at = 4 — 8). 

The next step is to carry out a similar investiga- 
tion for dynamical QCD. Many of the character- 
istics of DWF found here are sufficiently generic 
so that one would expect that they will also be 
present in QCD. If it turns out that QCD at the 
presently accessible lattice spacings, volumes and 
quark masses has similar restoration rates as the 
ones found here, then DWF will indeed provide a 
powerful fermion discretization method. 
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